18.997 Topics in Combinatorial Optimization February 12, 2004 

Lecture 4 

Lecturer: Michel X. Goemans Scribe: Constantine Caramanis 



This lecture covers: the Matching polytope, total dual integrality, and Hilbert bases. 

1 The Matching Polytope and Total Dual Integrality 

In this section we introduce the matching polytope as the convex hull of incidence vectors of match- 
ings. Next, we give a linear description of the matching polytope. and prove that the linear descrip- 
tion is correct by introducing the concept of Total Dual Integrality (TDI). 



Given a graph G = {V, E), and a matching M on G, we can identify M with its incidence vector: 



1 ee M 
otherwise 



We define the matching polytope V — V{G) to be the convex hull of these incidence vectors: 

P{G) = conv{x^ : M a matching of G.} 

We wish to obtain a linear description of V C IRI-^I. We must have > for e £ E. Also, every 
vertex can have at most one adjacent edge in any matching, and thus 

Thus our first attempt at a linear description is: 



pi _ f Xe>0 Ve e £■ 



x{S{v)) < 1 VvgV 

Consider the triangle, with edges labelled {1,2,3}. In this case, the matching polytope is 

P = conv{(0, 0,0), (1,0,0), (0,1,0), (0,0,1)}. 

The point (1/2,1/2,1/2) is in P^, i.e., it satisfies the constraints above, however it is not in the 
convex hull of the matching vectors. This example motivates the following family of constraints. 

Observe that for any matching M, any odd cardinality subset U can have at most {\U\ — l)/2 
edges. Thus we have the additional constraints 

x{E{U)) = J2 Xe< ^^^L^, UCV, \U\ odd, 

eeE{U) 

which we call the odd subgraph constraints. For the triangle, taking U to be the graph itself, we 
get the constraint xi + X2 + xs < 1. This constraint is violated by the point (1/2, 1/2, 1/2). As a 
second attempt at a linear description of the matching polytope we have: 

Xe>0 Ve e £ 

p2 = ^ xls{v)) < 1 Vw e F 

x{E{U)) < 1^ UCV, \U\ odd 
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Theorem 1 (Edmonds, 1965) The linear description is in fact the Matching polytope, i.e., 
V = P^. 



Proof: Edmonds gave an algorithmic proof: For any given weight vector, w, he argued that 
the optimization over gives the same maximum weighted matching as the actual solution. The 
algorithm uses shrinking and expanding of blossoms. It is in Chapter 27 of Schrijver's book. □ 
Rather than give Edmonds's algorithmic solution, we introduce the concept of Total Dual Integrality 
(TDI). This, in addition to Total Unimodularity, is further explored in Lecture 5. 
Recall the primal and dual standard formulation of an LP. 

(Primal (P)) | ^^^^ ■ '',^x<b\^{ ■ "^^^^^ " | ^^^^^ ^^^^ 
We define Total Dual Integrality as follows: 

Definition 1 (TDI) A polyhedron defined by a system of inequalities {x : Ax < b} (with A and b 
rational) is Total Dual Integral (TDI) if for any integral cost vector c to the primal problem, if the 
associated dual 

min : b^y 
s.t. : A^y = c 

y>0 

has a finite optimal value, then it also has an integral solution y* . 

We now state the main theorem about TDI polytopes. We save the proof for the next lecture. 

Theorem 2 (Edmonds-Giles, 1978) // the system {Ax < b} is TDI, and b is integral, then the 
polytope {Ax < b} is integral, i.e., all vertices are integral. 

Total Dual Integrality is a property of the representation of the polyhedron. To illustrate this point, 
consider the two dimensional polytope 




7' = conv{(0,3),(2,2),(0,0),(3,0)} 

This polytope may have many different representations. For example, 

r xi > 0, X2>0 ] 
r^< XI + 2x2 < 6 > 
1^ 2x2 + xi < 6 J 
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While the polytope is, evidently, integral, the system given is not TDI. For example, consider 
c= (1,1). The dual to 

min : d x 
s.t. : x^V 

is the linear program 

min : 6yi + 61/2 
s.t. : 2/1 + 2^2 > 1 

22/1 + 2/2 > 1 

2/1 , 2/2 > 

The dual optimal value is finite, however there are no integral optimal points. The dual tries to 
express (1, 1) as a positive linear combination of the vectors (1, 2) and (2, 1). This cannot be done 
as an integral linear combination. From this example, we see that in order for the dual to have an 
integral optimal solution, we must be able to express any integral vector in the dual cone C (as in 
the figure below) as an integral combination of some set of vectors. 




Definition 2 Given a rational (i.e. of the form {x : Ax < 0} with A rational) cone C, the minimal 
set of vectors such that any integral point in C can be expressed as an integral nonnegative linear 
combination of vectors in the set, is called a Hilbert Basis for the cone. 

For the cone defined by the vectors (1,2) and (2, 1), a Hilbert basis is given by the set of vectors 
H = {(1,2), (2, 1), (1, 1)}. We can get the additional vector (1, 1) by adding the redundant constraint 
xi + X2 < 4 in the primal. Therefore, the linear system 

Xi,X2>0 

x\ + 2x2 < 6 
2x\ + 3:2 < 6 

Xl + X2 < 4 

is TDI. 

Remark: Any rational cone has a finite Hilbert basis. This does not necessarily hold for all ir- 
rational cones. 

Theorem 3 For P a rational polyhedron, there exists a TDI system {Ax < b} such that P — {x : 
Ax < b}. Furthermore, the polytope P is integral if and only if we can take b to be integral. 

The primal form of the statement says that if P is a rational polyhedron such that maxjc^x : x £ P} 
is integral for every integral vector c, then the polyhedron is integral. 
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Suppose that for every rational supporting hyperplane (i.e., for every a, (3 such that a^x < (3 
for every x € P, and {a^x = /?} fl P ^ 0) there exists an integer point on the hyperplane. Then, 
the polytope is integral. 

This gives us an analog of Farkas's lemma: The system Ax = b has no integral solution if and 
only if A^y integral and b^y non-integral has a solution. 
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